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Abstract 

In this paper, we report theoretical calculations of the photo-induced absorption (PA) spectrum of 
naphthalene and anthracene, with the aim of understanding those excited states, which are invisible 
in the linear optical absorption. The excited state absorption spectra are computed from the l-B^ 
and the l-B^ states, and a detailed analysis of the many-body character of the states contributing 
to various peaks in the spectra is presented. The calculations are performed using the Pariser-Parr- 
Pople (PPP) Hamiltonian, along with the full configuration interaction (FCI) technique. The role 
of Coulomb parameters used in the PPP Hamiltonian is examined by considering standard Ohno 
parameters, as well as a screened set of parameters. The results of our calculations are extensively 
compared with the experimental data where available, and very good agreement has been obtained. 
Moreover, our calculations predict the presence of high intensity features which, to the best of our 
knowledge have not been explored earlier. We also present concrete predictions on the polarization 
properties of the PA spectrum, which can be verified in experiments performed on oriented samples. 

PACS numbers: 78.30.Jw, 78.20.Bh, 42.65.-k 
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I. INTRODUCTION 



Conjugated polymers have already replaced the traditional inorganic materials and 
emerged as popular choice for manufacturing opto-electronic devices, (author?)-'-'-'- The 
inter-molecular localized and the intra-molecular delocalized tc— electrons are credited for the 
interesting linear and nonlinear optical properties in these materials. The large nonlinear 
optical (NLO) response in such materials is favorable for various device applications such 
as optical switching, computing, and electro-opto modulation. (author?)-'-'- These useful 
applications have motivated scientists to study the excited states which are important for 
understanding the nonlinear optics of these materials. For example, in centrosymmetric 
molecules, states can be classified as even parity (gerade) and odd parity (ungerade) states. 
From the ground state, using linear optical spectroscopy, one can probe the low-lying odd 
parity excited states, (author?)- while the even parity states stay invisible. Photo-induced 
absorption (PA) technique gives an opportunity to investigate these even parity states, which 
contribute to the nonlinear optical processes in these materials. 

Among various conjugated polymers, oligo-acenes, due to their interesting electronic and 
optical properties, have come up as potential candidates for possible applications in the 
field of nonlinear optics. (author?)-'— Because of the alternant nature of these molecules, 
different states in addition to the point group symmetry (D 2h ) can also be classified ac- 
cording to their particle-hole parity (±). According to normal convention, we assign the 
ground state 1A 9 a negative (— ) particle-hole parity. As per the selection rules of D 2 h 
point group symmetry, for the ground-state linear optical absorption in oriented samples 
of oligo-acenes, two polarization channels are possible. Absorption of ^/-polarized photons 
leads the system to the B 2u manifold, while the x-polarized photons cause transitions to 
the B^ u states. (author?)- Same selection rules, therefore, also offer both these polarization 
possibilities for the excited state absorptions from the B 2u and B^ u states. Thus, from the 
B 2u type states, using x/^-polarized photons, B^ g / A~-type excited states can be reached, 
while from the B^ u states, using x/^-polarized photons, A~ / B^ g can be obtained (cf. Fig. 
(0D)). Here, the noteworthy point is the reverse symmetry-polarization relationship for the 
two types of manifolds. In the 70's, Bergman and Jortner(author?)— '— studied, these 
even parity excited states (A~ & B^ g ) using the two-photon absorption (TPA) spectroscopy. 
They measured the TPA spectra of naphthalene and anthracene in crystalline as well as 
in solution mode. Later, Dick and Hohlneicher measured the spectra over a wide range, 
between 3.59 eV to 6.20 eV for the naphthalene (author?)— and anthracene(author?)— 
molecules in solution. As far as the theory is concerned, the calculations have been 
performed using various methodologies like SDCI, (author?)— CASSCF, (author?)— '— 
MRMP,(author?)i£-iI LCAO-MO, (author?)^ exact PPP, (author ?)iaJtt CNDO/S-CI,-^ 
and CNDO-CI(author?)— on the molecular phases of these oligomers. To the best of our 
knowledge, none of the above mentioned calculations have reported the photo-induced ab- 
sorption spectrum of any of these materials. Moreover, most of the calculations are restricted 
to the low-order treatment of electron-correlation effects. Therefore, in this paper, we present 
a large scale correlated study of photo-induced absorption spectra of the first two members 
of acene family, namely, naphthalene and anthracene (cf Fig. ) - The PA spectra from 
the lB 2u and the states are computed, and a detailed analysis of the many-body 

character of the important excited states contributing to the spectra is performed. The 
calculations are done using the Pariser-Parr-Pople (PPP) model Hamiltonian along with 
full configuration interaction (FCI) technique. The standard, (author?) — as well as the 
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Figure 1: Schematic diagram of the essential states involved in the PA in oligoacenes, and their 
polarization characteristics. The arrows connecting two states imply optical absorption, with po- 
larization directions stated next to them. Location of states is not up to scale. 



screened parameters(author?)— are used in the PPP model Hamiltonian, and the results 
from both sets of calculations are compared with the experiments so as to analyze the role 
of Coulomb parameters. Low energy features of our calculated PA spectra match very well 
with the peaks observed in the two photon spectroscopy (TPS) based experiments.— dlMJA^. 
Moreover, our calculations predict the presence of intense absorption features in the higher 
energy range of the PA spectra, along with their polarization characteristics. Particularly, 
for the case of anthracene novel high intensity features well below the ionization threshold 
should be detectable in future experiments. Predicted polarizations properties of various 
features can also be tested in measurements performed on oriented samples. 

The remainder of this paper is organized as follows. In Sec. [Til we briefly review the 
theoretical methodology adopted in this work. In Sec. DTD we present and discuss our results 
for the photo-induced absorption from the l-B^u an d the 1B^ U excited states. A comparison 
is made with the experimental and the other available theoretical results. Finally, in Sec. 
IIVI we summarize our conclusions and provide prospects for future work. 



II. THEORY 

The schematic structures of naphthalene and anthracene are shown in FigEl Both the 
molecules are assumed to lie in the xy-p\ane with the conjugation direction taken to be along 
the x-axis. The carbon-carbon bond length has been fixed at 1.4 A, and all bond angles 
have been taken to be 120°. The reason of choosing this symmetric geometry, as against 
various other possibilities has already been discussed in our earlier paper. (author?)- Note 
that these structures can also be seen as two polyene chains of suitable lengths, coupled 
together along the ^-direction. 

The correlated calculations are performed using the PPP model Hamiltonian, which can 
be written as 

H = H Cl + H c , 2 + H Cl c 2 + H ee , (1) 
where H Cl and H C2 are the one-electron Hamiltonians for the carbon atoms located on the 
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(a) (b) 

Figure 2: Schematic drawings of (a) naphthalene and (b) anthracene, considered in this work 

upper and the lower polyene like chains, respectively. Hc 1 c 2 is the one-electron hopping 
between the two chains, and H ee depicts the electron-electron repulsion. The individual 
terms can now be written as, 

Hc x = —to ^2 Bk,k', (2a) 

(k,k>) 

Hc 2 = -to B v ( 2b ) 

and 

H Cl c 2 = ~t± ( 2c ) 

= U^2 n^nn + V itj {ni - l){n s - 1) (3) 

i i^j 

In the equation above, k, k' are carbon atoms on the upper polyene chain, /x, v are carbon 
atoms located on the lower polyene chain, while % and j represent all the atoms of the 
oligomer. Symbol (...) implies nearest neighbors, and B it j = J2a( c l a c j,°- + h.c). Matrix 
elements to, and t± depict one-electron hops. As far as the values of the hopping matrix 
elements are concerned, we took to = 2.4 eV for both intracell and intercell hopping, and 
t± = to consistent with the undimerized ground state for polyacene argued by Raghu et 
al. (author?)-^ 

The Coulomb interactions are parametrized according to the Ohno relation- 
ship, (author?)— 

^[///^(l + o.eim^// 2 , (4) 

where, K it j depicts the dielectric constant of the system which can simulate the effects of 
screening, U is the on-site repulsion term, and Rij is the distance in A between the ith 
carbon and the jth carbon. In the present work, we have performed calculations using "stan- 
dard parameters'™ with U = 11.13 eV and Kij = 1.0, as well as "screened parameters'— with 
U = 8.0 eV and = 2.0 (i ^ j) and K^i = 1.0. The screened parameters employed here 
were devised by Chandross and Mazumdar— with the aim of accounting for the interchain 
screening effects in phenylene based polymers. Our motive behind using these parameters 
is precisely the same. 

The starting point of the correlated calculations for the molecules are the restricted 
Hartree-Fock (HF) calculations, using the PPP Hamiltonian. The many-body effects beyond 
HF are computed using the full configuration interaction (FCI) method. Details of this CI- 
based many-body procedures have been presented in our earlier works.(author?)— 



From the CI calculations, we obtain the eigenfunctions and eigenvalues corresponding to the 
correlated ground and excited states of examined molecules. Using these eigenfunctions the 
dipole matrix elements amongst various excited states are computed. These dipole matrix 
elements, along with the energies of the excited states are, in turn, utilized to calculate 
various PA spectra. The number of excited states computed for each symmetry manifold is 
quite large, ranging typically from 70-100, allowing us to probe a broad energy range in the 
PA spectrum. 

III. RESULTS AND DISCUSSION 

In this section, we discuss our FCI results on the PA spectrum of naphthalene and an- 
thracene. Fig. ([3] and 2J) presents the PA spectra of naphthalene and anthracene, respec- 
tively. The spectra are computed using the standard (left panel) as well as the screened 
parameters (right panel) in the PPP Hamiltonian. For both 1B£ U and 1B^ U excited state 
absorption, we present the x- and y-polarized spectra separately, so as to facilitate compar- 
ison with experiments performed on oriented samples. The wave functions of the excited 
states contributing to various peaks in spectra are tabulated in tables (cf. Ref.—). Since, 
the calculations are performed using the FCI approach, they are exact within the model 
chosen and cannot be improved. Therefore, any discrepancy which these results may ex- 
hibit with respect to the experiments is a reflection of the limitations of the model or the 
parameters used, and not that of the correlation approach. At this point, we would like to 
note that Ramasesha and coworkers have also presented FCI results on naphthalene^ 1 ^ and 
anthracene.— 1 ^ However, for the case of naphthalene they only presented the energies of 
some of the low-lying states, but did not calculate transition dipoles. While, for anthracene 
they computed the transition dipoles among a few even and odd parity states, but the PA 
spectrum was not probed. Moreover, the energy range of the states computed by us extends 
higher than the calculations performed earlier, and therefore, can be used to interpret future 
experiments. 

In our earlier paper,- we showed that experiments as well as calculations suggest that for 
both the oligomers 1B£ U state is located at least 1.5 eV higher than 1B% U state. Therefore 
if in a PA experiment, the pump is tuned close to lB^ energy, using an appropriate probe 
beam one can investigate higher excited states which are dipole connected to Thus, 
we will measure excited state absorptions only from the 1B^ U state, without any inter- 
mixing from 1B^ U PA spectrum. Moreover, employing oriented samples, individual energy 
manifolds (A~/B^ g ) can be investigated by using x- or y-polarized probe beams. Similarly, 
one can measure the 1B^ U excited state absorption by tuning the pump energy to E(1B^ U ). 
Therefore, by choosing suitable pump energies and oriented samples, all the four components 
of the PA spectrum of each oligomer computed here, can be measured separately. We would 
like to mention here that 1B% U and 1B$ U PA spectra are plotted with respect to their energies 
obtained in linear optical spectra, reported in Ref. (-). In naphthalene, using standard 
parameters we calculated 1B} U state at 4.45 eV and 1B£ U state at 5.99 eV, while screened 
parameter values for these states were obtained to be 4.51 eV and 5.30 eV, respectively. In 
case of anthracene, standard parameters predicted EilB^y) =3.66 eV and E(1B^ U ) =5.34 
eV, on the other hand screened parameters calculated IB^ and 1B$ U states at 3.55 eV and 
4.64 eV, respectively. The energies reported in tables (cf. Ref.—) are the excitation energies 
of those states calculated with respect to the 1A~ ground state. 

We first present the lower end of the spectrum, which is that region where experimental 
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Figure 3: (color online) Photo-induced spectra of naphthalene computed using standard parameters 
(left panel) and screened parameters (right panel) in the PPP model Hamiltonian. The left side 
of both the spectrum depicts the PA from the IB^ state, while right side shows it from the IB^ U 
state. The x(y)-polarized transitions are shown using solid (dashed) curves. 

information about various excited states is available through two photon spectroscopy based 
experiments.— li^iiiiMiSi Although, the energies of these states are not so low in terms of their 
absolute values, but these states will appear in the low energy region of the PA spectrum. 
In absolute terms, the states for naphthalene and anthracene in this region are up to 6.5 eV 
and 5.2 eV, respectively, while in the PA spectra, they will appear in the range of 0-2 eV 
for naphthalene, and 0-1.6 eV for anthracene. Thereafter, we compare our results with the 
experimental data and calculations of other authors. Next, we present the discussion for the 
higher energy region of the spectrum which so far has not been probed, either experimentally, 
or theoretically. While discussing the high energy features we have to keep in mind that the 
ionization potentials of naphthalene and anthracene are 8.1 eV and 7.4 eV, respectively. 



A. Low energy region of the PA spectrum 

The first general aspect of the calculation is the relative location of the 2A~ state of 
these oligomers with respect to their 1B£ U state. Experimentally speaking both naphthalene 
and anthracene exhibit fluorescence, suggesting that 1B% U state is lower than 2A~ state. 
In our standard parameter calculations we obtain this ordering of states, while in screened 
parameter calculation 2A~ state is obtained below 1B£ U state. As a result, 2A~ state does 
not appear in the 1B% U PA spectrum, computed using the screened parameters, while in 
standard parameter spectrum it appears as the first (but very weak) feature in both the 
molecules (cf. left panel of Fig. (J3j) and Fig. (j4j)). After the 2A~ peak, the next important 
feature in the spectra corresponds to the l-B^ state, which is relatively more intense than 
the 2A~ feature. The intensity of the lB% g peak increases, while that of 2A~ decreases, with 
the increase in length of molecule, irrespective of the parameters used. Both features occur 
due to transition from the IB^ U state. Turning over to the 1B% U PA spectra, we notice that 
the low energy region gets very faint contributions from some of the states which contribute 
to the intermediate energy region of 1B% U PA spectra. 

Now we compare our results for naphthalene and anthracene with the experiments, and 
the works of other authors in a quantitative manner. For both the oligomers, this comparison 
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Figure 4: (color online) Photo-induced spectra of anthracene computed using standard parameters 
(left panel) and screened parameters (right panel) in the PPP model Hamiltonian. The left side 
of both the spectrum depicts the PA from the IB^ U state, while right side shows it from the 1B^ U 
state. The x/y-polarized transitions are shown using solid/dashed curves. 

is summarized in Table fl] For the case of naphthalene, our results obtained using the 
standard parameters are in very good agreement with the available experimental and other 
theoretical results, while screened parameters underestimate the experimental findings (cf. 
Tabled]). Our standard parameter calculations predict a very weak peak corresponding to 
2A~ state at 4.91 eV, while lB^ g , 3A g , and 4A~ states are obtained at 5.35 eV, 5.67 eV, and 
6.37 eV, respectively. Upon comparing our results to experiments, we find that the measured 
energy(author?)— '— '— of the 2A~ & 3A~ states (~5.50 eV & 6.05 eV) compares very well 
with our computed energy for 3A~ and 4A~ states (5.67 eV & 6.37 eV). Thus, the measured 
nA~ state consequently maps to (n + 1)A~ state of our work. We strongly believe that the 
2A~ state located around 4.91 eV (~39601 cm" 1 ) has been missed by the experimentalists due 
to its very weak intensity. This, in our opinion calls for further experimental investigations. 
While comparing theory and experiments for the anthracene, a similar mismatch in the 
labeling of the A~ states was pointed out by Ramasesha et al. (author?)— On the other 
hand, there is no ambiguity for the 1B% state, which is calculated at 5.35 eV by us and 
measured around 5.20 eV, using two photon spectroscopy (TPS). (author?)— '— '— Thus, 
the disagreement between our theory and experiments is no more than 5%. When this 
nA~-(n + 1)A~ mismatch is accounted for, we obtain very good agreement between our 
calculations and those of Hohmeicher and Dick, (author?)— who reported lB^ g , 2A g , and 
3A~ states at 5.31 eV, 5.45 eV, and 5.91 eV using SDCI/M approach, and at 5.37 eV, 5.70 
eV, and 6.12 eV using SDCI/P approach. Our results also agree well with the CNDO/S-CI 
calculations of Swiderek et al, (author?)— who computed lB^ g , 2A~ , and 3A~ states at 
5.37 eV, 5.61 eV, and 5.76 eV, respectively. On comparing our results with the MRMP 
and CASSCF calculations of Hashimoto et al., (author?)— we found that their CASSCF 
energies are much higher than ours, as well as compared to experimental energies, while 
MRMP results are comparable. Similarly, SCI results of Pariser(author?) 1B - predict much 
higher energies differences as compared to our values. 

In anthracene, the standard parameters predict 2A~ , lB% g , 3A~ , and 4A~ states at 3.89 
eV, 4.31 eV, 4.98 eV, and 5.14 eV, respectively, while the screened parameters values for 
the IB3 , 3v4~, and 4A~ states are 3.86 eV, 4.28 eV, and 4.61 eV, respectively. In order 
to compare our results with the two photon excitation spectrum measured by Dick and 
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Hohlneicher, (author?)— we found that they reported a very weak intensity peak around 
3.64 eV (~29300 cm" 1 ) and associated A~ symmetry to it. Although, they did not assign any 
state to this feature, yet from our computed PA spectrum and excitation energy it compares 
very favorably with calculated 2A~ state, a fact also noted by Ramasesha et al. (author?)— 
Therefore, following the lead of Ramasesha et al. (author?) 20 , we also compare experimental 
nA g states with our (n + l)A g states. This brings our standard parameters results in very 
good agreement with the experimental values of Dick and Hohlneicher, who reported lB^ , 
2A~, and 3A~ states at 4.44 eV, 4.71 eV, and 5.33 eV, respectively. (author?) 14 Prior to 
Dick and Hohlneicher(author?)~, Bergman and Jortner(author?)— also measured the two 
photon spectrum of anthracene molecule in solution, as well as of crystalline anthracene. As 
is obvious from Table [J that our standard parameter values are in good agreement with both 
their data. Since screened parameter results consistently underestimate the experimental 
values, we believe that they are not applicable for the case of anthracene as well. As far 
as comparison with other theoretical work is concerned, understandably our results are in 
excellent agreement with the work of Ramasesha et al, (author?)— who also used the FCI 
method in conjunction with the PPP model Hamiltonian. However, we have gone beyond 
the work of Ramasesha et al. (author?) 20 by also calculating the PA spectra of anthracene. 
On comparing our (n + l)A g states with the nA g states predicted by Kawashima et al.p- 
we obtain a very good agreement with the MRMP calculations. The calculations performed 
by other authors using SDCI/M, (author?)- 2 ^ SDCI/P, (author?) 22 CASSCF, (author?)^ 
and SCI, (author?) 18 approaches predict the states discussed at relatively larger energies 
than our results. 

B. Higher energy region of the PA spectrum 

This higher energy region, which includes the excited states in the range 6.6 eV to 8.1 
eV for naphthalene, and 5.3 eV to 7.4 eV for anthracene, has so far not been experimentally 
explored. In this region of the PA spectra highest intensity transitions correspond to the 
^-polarized photons, in agreement with the fact that the conjugation length of the systems 
is in the x-direction. Therefore, B 3g states lead to the most intense peaks in the 1B£ U PA 
spectra, while A~ states contribute to the main intensity of the IB^ U PA spectra. However, 
when these two x-polarized transitions are compared to each other, B 3g peaks are found to 
be more intense than A~ peaks. Thus, B 3g states which exhibit strong dipole coupling to 
1-Bjj may end up playing the same role in the nonlinear properties of these materials which 
the mA~ state plays for trans-polyacetylene.— ^ 

For naphthalene, with standard parameters 2B 3g state (6.58 eV) leads to low intensity 
features (c/. Fig. [3j), while the 3B 3g state (8.1 eV) close to the ionization threshold leads 
to a very intense peak in IB^ U PA spectrum. Similarly, 6A g (8.06 eV) also located close 
to the ionization threshold contributes a very intense peak to the IB^ U PA spectrum and 
a weak one to the 1B% U spectrum. Because of its proximity to the ionization threshold it 
may be difficult to detect the 6A~ and 3B 3g states, but the 2B 3g state should certainly be 
detectable. 

The PA spectrum of anthracene computed using the standard parameters is much richer 
compared to that of naphthalene, and offers many more detectable peaks. For example, in 
the B^g manifold, states 2B 3 ~ g (6.37 eV), 3B 3g (6.68 eV), and AB 3g (6.98 eV), are well below 
the ionization threshold of 7.4 eV and lead to significant peaks either in the lB^ or lB 3u 
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PA spectra. Particularly noteworthy is the strong intensity of the 3B 3g peak in the PA 
spectrum. 5B 3g state at 7.39 eV also leads to an intense peak but may be difficult to detect 
because of its closeness to the ionization threshold. As far as the A~ manifold is concerned, 
bA~ (5.77 eV), 6A~ (6.45 eV), and 7A~ (7.16 eV) states contribute significant intensity to 
the PA spectrum and have not been explored so far. Reasonablly strong intensities, coupled 
with the fact that they are away from the ionization threshold makes them strong candidates 
for experimental investigation. The 8A g (7.3 eV) peak which is also quite intense is close to 
the ionization threshold and hence may be difficult to detect. 

As far as the screened parameters based results are concerened, for both the oligomers 
more features are visible in the same energy range, because of the fact that the spectrum is 
red shifted as compared to the standard parameters. This trend is in agreement with our 
earlier work on linear optical absorption in these materials.- For example, in naphthalene 
states all the way up to 45^ and 7A g , while for anthracene states up to 6B 3g and 1L4.J, 
are below the ionization threshold of individual oligomers, and therefore, in principle are de- 
tectable. Thus, the experimental invesitgation of these high energy features can throw light 
on the nature of Coulomb interactions in these materials. Smaller number of peaks below 
the ionization threshold will indicate that the standard parameters describe the Coulomb 
interactions in these systems, while the larger number of features will suggest the validity 
of the screened parameters. 

Regarding wave functions, one can notice from the tables (cf. Ref.— ), that single as well 
as double excitations contribute equally to various excitated states visible in these spectra. 
In particular, the main configuration contributing to the many-particle wave function of the 
2A~ state corresponds to a double excitation described by HOMO (iJ)^LUMO (L); HOMO 
(i/)^LUMO (L) (\H -> L;H -> L ) ), while single excitation \H -> L + 2 ) ± c.c. is found 
as the major contributor to the many-particle wave function of the 1-RT state. The most 
intense B 3g feature for both naphthalene and anthracene computed using the standard as 
well as the screened parameter occurs mainly due to \H — > L; H — > L + 1 ) ±c.c. excitation. 
For the case of A~ peaks, maximum intensity can be attributed to states which for both the 
oligomers, and irrespective of parameters used, are described by \H — > L; H — 1 — ► L + 1 ) 
configuration. 

IV. CONCLUSIONS 

In this paper, we presented FCI method based large scale correlated calculations of photo- 
induced absorption spectrum of naphthalene and anthracene from their low-lying one-photon 
states, 1B£ U and 1B£ U . The aim behind these calculations was to understand the nature of 
important two photon states of these materials and their coupling to one photon states. The 
spectra were computed over a broad energy range, so far unexplored in earlier calculations 
as well as experiments. The main conclusions of this work are as follows: 

1. Our most reliable results were obtained using standard Coulomb parameters in the 
Hamiltonian, leading to excellent agreement on the energies of low lying two photon 
states, available from two photon spectroscopy (TPS) experiments. 

2. In particular, the 2A~ state was calculated to be higher than the 1B£ U state, predicting 
these materials to exhibit fluorescence, in agreement with experimental data. However, 
two photon spectroscopy based experiments have not been able to locate the 2A~ state, 
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Table I: Comparison of results of our calculations performed with the standard (Std.) parameters 
and the screened (Scd.) parameters with other experimental and theoretical results for the most 
important low-lying even parity states. The energies of various states are given with respect to the 
ground state (IA~). In naphthalene and anthracene, using the standard parameters IB^ U state was 
obtained at 4.45 eV and 3.66 eV, respectively- The screened parameters predicted 1-B^ at 4.51 eV 
for naphthalene, and 3.55 eV for anthracene.— 

Excitation energy (eV) 

State Present work Experimental Other theoretical 

Std. Scd. 





para 


para. 




Naphthalene (CioHg) 


2A- 


4.91 










5.35 


4.77 


5.22, (author?)^ 5.21, (author?)^- 5.27(author?)^ 


5.31 a , (author?) 22 5.37 b , (author?) 




5.67 


5.01 


5.52, (author?)^ 5.45,(author?)JJ. 5.50(author?)^ 


5.45 a , (author?)^ 5.70 b , (author?) 


4A- 


6.37 


5.62 


6.05(author?)i2 


5.91 a , (author?) 22 6.12 b , (author?) 


Anthracene (C14H10) 


2A- 


3.89 




3.64 e (author?)^ 






4.31 


3.86 


4.44 e , (author?)^ 4.53 f , (author?)^ 4.66 e , (author?)^ 


4.31, (author?)^ 4.57 a , (author?) 2 


3A~ 


4.98 


4.28 


4.71 e , (author?)^ 4.81 f ,(author?) 12 4.96 e , (author?)^ 


4.96, (author?)^ 4.29 a , (author?) 5 


4A- 


5.14 


4.61 


5.33 e (author?)^ 


5.12, (author?)^ 5.18 a , (author?) 1 


a SDCI/M ca 


culations d CASSCF calculations 



b SDCI/P calculations e solution spectrum 

c MRMP calculations Crystalline spectrum 

possibly due to its very low intensity. Therefore, we encourage careful TPS or PA based 
measurements of the 2A~ state in these materials to verify our predictions. 

3. We make specific predictions about the higher energy region of the PA spectra, where 
we found intense features corresponding both to the A~ and B^ g type states. These 
states may also be visible in the TPS experiments in the higher energy region, which so 
far has not been explored. We further predict that the relative intensities of the peaks 
corresponding to x-polarized transitions is more than the y-polarized ones, which can 
be probed in experiments performed with the oriented samples. 

4. We believe that it is feasible to resolve various components of the PA spectra computed 
here, by tuning the pump beam to different energies {EiAB^y) or E{1B^ U ))^ and by 
using oriented samples. 

Because of all the reasons mentioned above it is of interest to probe the PA spectra of these 
well established materials; therefore, we urge strong experimental efforts in this direction. 

The calculation of TPA spectrum will additionally verify our results, as the peaks appear- 
ing in the PA spectra should also appear in the TPA spectra of the considered oligomers. 
Moreover, some of the important features of the PA spectra may also appear in the third 
harmonic generation (THG) spectra, giving more insight into the states governing the non- 
linear optics of these materials. Therefore, in near future we will present the results of 
calculations of the TPA and THG spectra of oligo-acenes, which are currently underway in 
our group. 
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SUPPLEMENTARY INFORMATION 



Here we present the tables summarizing the results of our FCI calculations for excited 
state absorption in naphthalene and anthracene. The data presented in the tables includes 
important configurations contributing to the many-body wave functions of various excited 
states, their excitation energies, and transition dipoles connecting them to the lB^ state 
and lBt, state. 



A. Naphthalene 



Table II: A~-type excited states contributing to the photo- 
induced absorption spectrum of naphthalene, corresponding 
to transition from IB^ (at 4.45 eV) and IB^ U state (at 5.99 
eV) due to the absorption of y-polarized and x-polarized pho- 
tons, respectively, computed using the FCI method coupled 
with the standard parameters in the PPP model Hamiltonian. 
The table includes many particle wave functions, excitation 
energies, and dipole matrix elements of various states with 
respect to IB^ and IB^ U states. '+c.c.' indicates that the 
coefficient of charge conjugate of a given configuration has 
the same sign, while '— c.c.' implies that the two coefficients 
have opposite signs. 



State 


E (eV) 


Transition Dipole (A) 


Wave Functions 






y- component 


x- component 


2A- 


4.91 


0.162 




\H -2 -»• L + l) +c.c.(0.4120) 










\H -» L;H -> L) (0.4047) 










\H -> L + 3) +c.c.(0.3555) 


3A- 


5.67 


0.147 




\H -» L; H - 1 -> L + 1 } (0.4062) 










\H - 1 -»• L + 2) - c.c.(0.3918) 










\H -► L + 3) +c.c.(0.3253) 










\H -» L;H -> L) (0.2196) 




6.37 


0.697 


0.242 


\H -» L;H -> L) (0.6117) 










\H -» L;H - 2 -> L + 2 } (0.2732) 










\H -» L; # - 1 -> L + 1 } (0.2556) 










|Jf -► L + 2;H -» L + 2) + c.c.(0.2537) 


6A- 


8.06 


0.300 


0.734 


|Jf -» L; # - 1 -> L + 1 } (0.4981) 










\H -> L;H -2-> L + 2) (0.3886) 










\H ^ L + 2;H ^ L + 2) + c.c.(0.2912) 










|#-2^L + 4) + c.c.(0.2143) 
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Table III: excited states contributing to the photo- 

induced absorption spectrum of naphthalene, corresponding 
to transition from l-B^ (at 4.45 eV) and 1B^ U state (at 5.99 
eV) due to the absorption of x-polarized and y-polarized pho- 
tons, respectively, computed using the FCI method coupled 
with the standard parameters in the PPP model Hamiltonian. 
The table includes many particle wave functions, excitation 
energies, and dipole matrix elements of various states with 
respect to lB^ and IB^ U states. '+c.c.' indicates that the 
coefficient of charge conjugate of a given configuration has 
the same sign, while '— c.c.' implies that the two coefficients 
have opposite signs. 



State 


E (eV) 


Transition 


Dipole (A) 


Wave Functions 


x- component 


y- component 


IB* 


5.35 


0.672 




\H ^L + 2) +c.c.(0.5469) 










\H -► L; H -» L + 1 ) + c.c.(0.2841) 


2^ 


6.58 


0.154 


0.219 


\H - 1 -► L + 3 } + c.c.(0.4828) 










\H -1^ L + l-H ^ L + l) + c.c.(0.3665) 




8.12 


1.271 


0.255 


\H -> L; H -» L + 1 ) + c.c.(0.4081) 










\H ^L + 2) +c.c.(0.2349) 










\H -► L;H -» L + 4) + c.c.(0.2335) 










\H -?> -> L;H -> L + 2) + c.c.(0.2256) 










\H -2^L) + c.c.(0.2128) 



Table IV: A~-type excited states contributing to the photo- 
induced absorption spectrum of naphthalene, corresponding 
to transition from IB^ (at 4.51 eV) and 1B% U state (at 5.30 
eV) due to the absorption of y-polarized and x-polarized pho- 
tons, respectively, computed using the FCI method coupled 
with the screened parameters in the PPP model Hamiltonian. 
The table includes many particle wave functions, excitation 
energies, and dipole matrix elements of various states with 
respect to l-B^u an d l-SjL states. '+c.c.' indicates that the 
coefficient of charge conjugate of a given configuration has 
the same sign, while '—c.c' implies that the two coefficients 
have opposite signs. 



State 


E (eV) 


Transition 


Dipole (A) 


Wave Functions 


y- component 


x-component 




5.01 


0.355 




\H -» L; H - 1 -» L + 1 ) (0.4022) 










\H -► L + 3) -c.c.(0.3582) 










\H- 1 -^L + 2) +c.c.(0.3396) 




5.62 


0.856 


0.612 


\H -> L;H -» L) (0.5903) 


Continued on next page 
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Table IV - continued from previous page 



State 


E (eV) 


Transition Dipole (A) 


Wave Functions 


y- component 


x-component 










\H -» L;H - 2 -» L + 2) (0.2873) 










\H- 1 -^L + 2) +c.c.(0.2594) 










|if -» L + 2-H -> L + 2) - c.c.(0.2494) 


5A~ 


7.25 


0.336 


1.331 


|fl'->L;fl'-l-»L+l) (0.5808) 










|fr->L+l;fl"^L + l> - c.c.(0.3109) 


6A- 




- 


0.872 


|fl'->L;fl'-l-»L+l) (0.3260) 










|if -2^L + 4) + c.c.(0.3156) 










\H -> L;H -2 -> L + 2) (0.2960) 










|if -» L + 2;H -> L + 2) - c.c.(0.2204) 


7A- 


8.13 




0.657 


|if-4^L + 2) + c.c.(0.2293) 










|fl"->L;fT-l-»L+l) (0.2262) 


8A- 


8.48 


0.978 


0.817 


|fT-l->L + l;fr-l-»L+l) (0.3318) 










|if-2^L + 2;if-l^L + l) (0.2959) 










|if -» L; if - 1 -» L + 1 ) (0.4022) 










|ff ->■ L + 3) - c.c.(0.2825) 



Table V: 5^-type excited states contributing to the photo- 
induced absorption spectrum of naphthalene, corresponding 
to transition from \B^ U (at 4.51 eV) and 1B^ U (at 5.30 eV) 
states due to the absorption of x-polarized and y-polarized 
photons, respectively, computed using the FCI method cou- 
pled with the screened parameters in the PPP model Hamil- 
tonian. The table includes many particle wave functions, ex- 
citation energies, and dipole matrix elements of various states 
with respect to l-B^ and IB^ U states. '+c.c.' indicates that 
the coefficient of charge conjugate of a given configuration has 
the same sign, while '— c.c.' implies that the two coefficients 
have opposite signs. 



State 


E (eV) 


Transition 


Dipole (A) 


Wave Functions 


x- component 


y- component 


IB* 


4.77 


1.119 




\H ^L + 2) + c.c.(0.5519) 










|if -» L; H -» L + 1 ) - c.c.(0.2247) 




5.75 




0.263 


|ff - 1 -> L + 3) +c.c.(0.4766) 










\H-1^>L+1;H->L + 1) + c.c.(0.3324) 




7.09 


1.557 


0.259 


|ff -» L; H -» L + 1 ) - c.c.(0.4407) 










|ff -► L + 2} - c.c.(0.2124) 










|ff^L + 2;ff^L + 3} + c.c.(0.2044) 




7.94 


1.194 


0.159 


|if -► L + 2;H -» L + 3) + c.c.(0.2738) 










\H —> L; H —> L + 4) - c.c.(0.2424) 










|if -» L; if -» L + 1 ) + c.c.(0.2418) 


Continued on next page 
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Table V — continued from previous page 



State 


E (eV) 


Transition Dipole (A) 


Wave Functions 


x- component 


y- component 










\H + +c.c.(0.2089) 




8.39 


0.267 


0.194 


\H-4 -> L + 3) + c.c.(0.3131) 










\H -► L + 3;H - 1 ^ L + 3} - c.c.(0.2058) 
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B. Anthracene 



Table VI: ^4~-type excited states contributing to the photo- 
induced absorption spectrum of anthracene, corresponding to 
transition from l-B^u ( a ^ 3.66 eV) and 1-B^ state (at 5.34 
eV) due to the absorption of y-polarized and x-polarized pho- 
tons, respectively, computed using the FCI method coupled 
with the standard parameters in the PPP model Hamiltonian. 
The table includes many particle wave functions, excitation 
energies, and dipole matrix elements of various states with 
respect to l-B^ and 1-B^ states. '+c.c.' indicates that the 
coefficient of charge conjugate of a given configuration has 
the same sign, while '— c.c.' implies that the two coefficients 
have opposite signs. 



State 


E (eV) 


Transition 


Dipole (A) 


Wave Functions 


y- component 


x-component 


2A~ 


3.89 


0.150 


- 


|ff -» L;H -> L) (0.5639) 










|if -» L + 3) +c.c.(0.3372) 










\H - 1 -> L + 2 } + c.c.(0.2476) 


3A~ 


4.98 


0.635 


- 


|ff -» L + 3) - c.c.(0.3502) 










|if -> L; if -1 -> L+l) (0.3331) 










|if -> L;if -> L) (0.3015) 










\H ^ L;H - 2 ^ L + 2) (0.2562) 


4A~ 


5.14 


0.434 




|fl"-l-»L + 2> + c.c.(0.4913) 










\H ^ L;H ^ L) (0.3685) 










|if ->■ L; if - 1 -> L + 1 } (0.2605) 


5A- 


5.77 


0.147 


0.201 


|ff ->■ L;H -» L + 4 ) - c.c.(0.3709) 










|if -» L + 6) + c.c.(0.3412) 










\H ^ L;H - 2 ^ L + 2) (0.3095) 


6A- 


6.45 


0.304 


0.233 


\H ^ L;H - 2 ^ L + 2) (0.3396) 










|ff -»• L + 2;H -» L + 2) - c.c.(0.3212) 










|if -»• L; if - 1 -> L + 1 } (0.2204) 


7A~ 


7.16 




0.627 


|fl'->L+l;fl'-fL + l> + c.c.(4153) 










|ff -► L; if - 1 -> L + 1 } (0.4150) 


8A~ 


7.30 


0.155 


0.691 


|ff -► L; if - 1 -> L + 1 } (0.3325) 










|if-4^L + 3} + c.c.(0.2810) 










|fr->L + 2;fl"-l-»L + 3) + c.c.(2137) 










|ff^L;if-4^L + 4) (0.2027) 


10A- 


7.92 


0.348 


0.403 


|ff -> L; if - 1 -> L + 1 } (0.2848) 










|ff -► L;ff - 2 -»• L + 2) (0.2178) 
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Table VII: ff^-type excited states contributing to the photo- 
induced absorption spectrum of anthracene, corresponding to 
transition from IB^ U (at 3.66 eV) and IB^ U state (at 5.34 
eV) due to the absorption of x-polarized and y-polarized pho- 
tons, respectively, computed using the FCI method coupled 
with the standard parameters in the PPP model Hamiltonian. 
The table includes many particle wave functions, excitation 
energies, and dipole matrix elements of various states with 
respect to lB^ and 1B^ U states. '+c.c.' indicates that the 
coefficient of charge conjugate of a given configuration has 
the same sign, while '— c.c.' implies that the two coefficients 
have opposite signs. 



State 


E (eV) 


Transition 


Dipole (A) 


Wave Functions 


x-component 


y-component 


IB* 


4.31 


1.037 




|ff -► L + 2) +c.c.(0.5226) 










|ff -» L; ff -► L + l) + c.c.(0.2808) 




6.37 




0.278 


|ff -► L + 2) - c.c.(0.6268) 




6.68 


1.168 


0.126 


|ff-l^L + 3) - c.c.(0.3000) 










|ff -» L; H -> L + 1 } + c.c.(0.2970) 




6.98 




0.155 


|ff -4 -» L + 2) +c.c.(0.3138) 










|ff->L + l;ff-»L + 4) + c.c.(0.2492) 




7.39 


1.253 


0.245 


|ff-2^L + 2;ff^L + l) + c.c.(0.2496) 










|ff ^ L; ff ^ L + 1 ) - c.c.(0.2449) 










|ff ^L;ff-l^L + 4) +c.c.(0.2101) 










|ff -4 -» L + 2) +c.c.(0.2009) 




7.77 


1.313 


0.329 


|ff -» L; ff ->■ L + 1 ) - c.c.(0.2409) 










|ff^L + 2;ff^L + 3) + c.c.(0.2293) 



Table VIII: ^4~-type excited states contributing to the photo- 
induced absorption spectrum of anthracene, corresponding to 
transition from lff^ (at 3.55 eV) and 1B^ U state (at 4.64 
eV) due to the absorption of y-polarized and x-polarized pho- 
tons, respectively, computed using the FCI method coupled 
with the screened parameters in the PPP model Hamiltonian. 
The table includes many particle wave functions, excitation 
energies, and dipole matrix elements of various states with 
respect to lff^ and 1B^ U states. '+c.c.' indicates that the 
coefficient of charge conjugate of a given configuration has 
the same sign, while '—c.c' implies that the two coefficients 
have opposite signs. 



State 


E(eV) 


Transition 


Dipole (A) 


Wave Functions 


y- component 


x- component 


Continued on next page 
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Table VIII - continued from previous page 



State 


E (eV) 


Transition Dipole (A) 


Wave Functions 


y- component 


x- component 


3,4- 


4.27 


0.782 


- 


\H -► L + 3) +c.c.(0.3743) 










\H -» L; H - 1 -> L + 1 ) (0.3083) 










\H -> L;H -> L) (0.2538) 










\H -» L;H - 2 -> L + 2 ) (0.2284) 


4,4" 


4.62 


0.461 


- 


|Jf -»• L + 3) +c.c.(0.4822) 










\H ^ L;H ^ L) (0.3492) 










|fl'->L;fl'-l-fL + l) (0.2122) 


5A~ 


5.28 


0.139 


- 


|JJ -» L; # -» L + 4) - c.c.(0.3492) 










|#^L + 6) - c.c. (0.3427) 










|if -» L; # - 2 -> L + 2) (0.2860) 


6A- 


5.72 


0.536 


0.453 


|if -» L; # - 2 -> L + 2) (0.3011) 










|#^L + 2;#^L + 2) + c.c.(0.2867) 


7A~ 


6.45 


0.178 


1.767 


|fl"->L;fl"-l-fL + l) (0.5117) 










|L 7 -3^L + 4) -c.c.(0.2266) 










|fr->L+l;fl"-fL + l> + c.c.(0.2084) 


9A~ 


6.62 


- 


1.038 


|#-4^L + 3) + c.c.(0.2661) 










|#^L;#-1^L + 1) (0.2513) 


WA- 


7.11 


0.373 


0.403 


|if -» L-H - 2 -> L + 2) (0.2783) 










|#^L;#-1^L + 1) (0.2491) 










|F-2^L + 2;#^L + 4) + c.c.(0.2076) 


11A- 


7.34 


0.834 


1.437 


|iT -»• L + 3) +c.c.(0.2264) 










|#-2^L + 5) + c.c.(0.2144) 










|#^L;iI-l^L + l) (0.2128) 










\H -> L;H -> L) (0.2086) 


12A- 


7.63 


0.206 




\H ^ L + 1;H -3^ L + 2) + c.c.(0.2416) 










|ff-l->L+l;fl"-2-»L + 2) (0.2401) 
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Table IX: B^-type excited states contributing to the photo- 
induced absorption spectrum of anthracene, corresponding to 
transition from IB^ U (at 3.55 eV) and IB% U (at 4.64 eV) states 
due to the absorption of x-polarized and y-polarized photons, 
respectively, computed using the FCI method coupled with 
the screened parameters in the PPP model Hamiltonian. The 
table includes many particle wave functions, excitation ener- 
gies, and dipole matrix elements of various states with respect 
to l-B^u an d states. '+c.c.' indicates that the coefficient 
of charge conjugate of a given configuration has the same sign, 
while '— c.c.' implies that the two coefficients have opposite 
signs. 



State 


E (eV) 


Transition 


Dipole (A) 


Wave Functions 


x- component 


y-component 


IB* 


3.86 


1.593 




\H -» L + 2) +c.c.(0.5375) 










\H -» L; H -► L + 1 } - c.c.(0.2162) 


2^ 


5.52 


0.971 


0.373 


\H L + 2,) + c.c.(0.3189) 










\H L + l-H -> L + l) + c.c.(0.2503) 










\H -» L; H -► L + 1 } - c.c.(0.2053) 




5.81 


1.612 


0.140 


\H -» L; if ->■ L + 1 } - c.c.(0.3673) 










|if- 1 -» L + 3) + c.c.(0.3316) 




6.21 


0.613 


0.153 


|ii-4^L + 2) + c.c.(0.4228) 




6.43 


1.271 


0.174 


\H -» L; if -► L + 5) + c.c.(0.2224) 










|ii-2^f + 2;if^L + l) + c.c.(0.2040) 




6.75 


1.050 


0.220 


|if^L + 2;ii^L + 3} + c.c.(0.2451) 
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